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Abstract. We consider the nonequilibrium time evolution of the transverse magnetization in the critical 
Ising and XX quantum chains. For some inhomogeneously magnetized initial states we derive analytically 
the transverse magnetization profiles and show that they evolve into scaling forms in the long-time limit. 
In particular it is seen that the Ising chain exhibits some similarities with the conserved dynamics XX 
chain. That is, after a transient regime, the total residual magnetization in the transverse direction is also 
conserved in the Ising case. A class of general initial states is also considered. 
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PACS. 75.40.Gb Dynamic properties (dynamic susceptibility, spin waves, spin diffusion, dynamic scaling, 
etc.) 

05.70.Ln Nonequilibrium and irreversible thermodynamics 
05.30-d Quantum statistical mechanics 
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1 Introduction 

Nonequilibrium properties of quantum systems have at- 
tracted a lot of interest since they have natural dynam- 
ics in contrast to classical ones and since classical effects 
are actually quantal. There are several ways to consider 
nonequilibrium quantum systems. One is to couple the 
quantum mechanical system to a heat bath which can be 
itself described quantum mechanically [Q . In this case, a 
part of the whole system is isolated and called 'the system' 
while the rest is supposed to describe a certain environ- 
nement with which the system interacts and dissipates 
through. Another route is to impose a current on the sys- 
tem and investigate the steady states Still an other 
possibility is simply to investigate the relaxation of an 
initial state, in which the system has been prepared, and 
evaluate expectation values of observables at later times. 
This was done recently on the XX-quantum chain 
with a step-like magnetization initial state. More recently 
in ref. Q , the relaxation of spatially inhomogeneous ini- 
tial states has been treated for several variants of the XY 
quantum model. Relaxation phenomena at zero tempera- 
ture with homogeneous initial state has been considered 
in ref. for the XX and Ising chains in a transverse 

field in the context of aging. 

In this work we study the nonequilibrium profiles of the 
critical Ising and XX quantum chains. We suppose that 
at the initial time t = the system is prepared in a given 
state The time evolution of the system is entirely gov- 
erned by the Schrodinger equation and is formally given 

by 

\^{t)) = expi^iHtm (1) 



since the systems under consideration are closed. The 
basic quantity we calculate is the expectation value of 
the local transverse magnetization at time t, m{l,t) = 
{ip{t)\(Tf\'I'{t)). We first consider two different initial states 
a kink in the z direction, | . . . TTTiii • • •)> which was 
already considered in ref. |^ for the quantum XX chain, 
and a droplet configuration L. . . It | .. Itt . . .) for both 
Ising and XX chains. In ref. |Q the cumulated magnetiza- 
tion was calculated for the XX chain but the profile itself 
was not considered. We give finally the general expression 
for the relaxation of the transverse magnetization in terms 
of a convolution product in the continuum limit. The ker- 
nels of both XX and Ising chains are readly expressed in 
both direct and Fourier space. 



2 Basic Quantities 

The one dimensional Ising and XX Hamiltonians with L 
sites and open boundary conditions are given by the same 
one-parameter anisotropic XY Hamiltonian: 



H = - 



1 

-Y 

k=l 



1 



1- K. 



'k"k+l 



fc=i 
(2) 

where the anisotropy parameter k = 1 corresponds to 
the Ising case with a Z2 symmetry and k = de- 
scribes the XX -model which has U{1) symmetry. The 
Hamiltonian (||) is diagonalizable through of a Jordan- 
Wigner transformation, followed by a canonical transfor- 
mation 10]. In terms of the Clifford operators {-T/}, the 
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Jordan- Wigner transformation is expressed as 



r„^--(-i)"-Mn"=>5 



(3) 



The generators {F^} satisfy 



in,l^l,...,L; i,j = l,2), (4) 



where we have introduced a pseudoscalar product defined 
as {C\D) = ^{C,D} with {., .} the anticommutator. The 
original spin variables are obtained in terms of the Fs by 
inverting the previous relations. One obtains cr^cr^+i — 



that (0) i 



rVn-y 



1; 

is written as 



-ir^+iF^ and 



-ir^rl so 



H = 



4 ^ 

k 



^rl{iay)rk+i+Kryrk+i-hrliayrk , (5) 



where = {r^,r^), the hermitian conjuguate of Fk, 
is a 2-components spinor and and cr*' are the Pauli 
matrices. Introducing the 2L-components Clifford opera- 
tor r^' = (r/, ^2^ . . . , rl), we arrive a.t H = (l/4)rtTr 
with T''' — T. The diagonalisation is then performed by 
the introduction of the diagonal Clifford generators 7^ = 
(7^,7^) related to the old one by — J2q4'qi^)lq ^^'^ 
rf = i{}q{l)^'^ with real (j) and tp components. Introduc- 
ing the Fermi operators r]q = 1/2(7^-1-17^) and the adjoint 
finally one arrives at the usual free fermionic Hamilto- 
nian H = £qri^riq + Eq. The excitations energies Eg and 
the transformation coefficients (j)q{l),^q{l) are solution of 
the 2L X 2L eigenvalue system TV^ ~ ^q^qi with com- 
ponents V^{k) = {(j)q{k),—itpq{k)). The eigenvectors sat- 
isfy the orthogonality relations J2q4'q{''-)4'qU) = ^ij and 

The time evolution of the spin operators are easily ex- 
pressed in terms of the time dependence of the Clifford 
generators F (see Appendix) . The basic time evolution of 
the diagonal operators is jq{t) = e*^*7ge~'^* = TZ{eqt)jq 
where TZ{0) is a rotation of angle 9. In matrix form we 
have 



cos eqt sin eqt 
— sinegi cos eqt 



(6) 



Using this relations, we can express the time dependence 
of the Fs through an expansion onto the basis {FJ;.}: 



JHt 



Fie 



-iHt 



k.i 



inim)) n (7) 



with components 

(Al|r/(t)) =^0g(fc)0g(Z)COSegt 

{Fl\F^{t)) - {F^\Fli-t)) =-Y,'Pq{k)i,q{l)sineqt 
{rl\Ff(t)) ^J2Mk)Ml) cos eqt . (8) 



For the Ising chain, at the critical point h = 1, 
the basic contractions are obtain in a closed form. For 
open boundary conditions, the excitation energies eq = 
2sin((7/2) and eigenvectors </> and ip are 



Ml) - (-1)' 



V2L+1 

Ml) = (-1)'+^^ 



V2L+1 



cos q{l - 1/2) 



sin ql 



(9) 



with q — (2p + 1)tt/{2L + 1). In the thermodynamic limit 
L oo, the contractions are then expressed in terms of 
Bessel functions J„(z) of integer order as ||ll[: 

(F^lFlit)) {F^\Fnt)) = (-1)'^+'J2(,_,)(20 

inirnt)) = -(-i)'^+'+v2(z-,)+i(2<) . (10) 

For the XX-chain in a similar way one obtains 

cos{ht) ; I — k = 2p 



{F^\Fnt)) = {^y-'Jl-k{t) 



—i sm{ht) ; I — k = 2p+ 1 



<Lnd(r^|r,'(t)) = (-r^!|r,Hf)). 

In what follows we are interested in the time relaxation 
of the transverse magnetization. In the Heisenberg picture, 
we have, using the expansion (0), erf (i) — —iF^{t)Fl{t) = 
-^Effk2in\\rnt)){nl\rHt))Fl\Fll The expectation 
value in the z-direction \'P') state is then 

m(M)=E[(A'l^'W)(A'l^''W) 

k 

- {F^\F^{t)){F^\Fnt))]{^Wm (12) 

since only the terms {W\F^F^\^) are non-vanishing in the 
state I If'). 



3 Kink-like initial state 

We consider first the initial state with a kink located at 
the origin 

(13) 



where | f, j) are the eigenstates of the Pauli matrix, 
"■^1 Tj i) = i| T: i)j and we take the thermodynamic limit 
N ^ oo. We present here only the results for the criti- 
cal {h — 1) Ising quantum chain since the XX chain was 
already considered in ref. Using the previous formu- 
las (p^ and (10) for the contractions, a straightforward 



calculation leads to 



m{l,t)^- J2 
p=i-i 



1 - 



Jim + j^r, m 



(14) 
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Now the analysis proceeds along the same lines as in 
ref. HI- We introduce the discrete derivative 



^n(«) = -t[m{n +l,t)- m{n, t)]„/t=„ 

= -2^[{l~v^)JU2t) + j'J{2t)] 



(15) 



For obvious symmetry reasons, we will consider only the 
part u > since we have m{—n,t) = —m{n,t). Due to 
the different asymptotic properties of the Bessel functions 
one has to distinguish between the cases v > 1 and v < 1. 
For V > 1, which means that n > t, we are in the acausal 
region (outside the light-cone) since the excitations, trav- 
eling with velocity one |10 , have no time to propagate 



from the initial position of the kink to the site n. Then 
it is the local environnement which completly governs the 
behavior of the magnetization. That is, the magnetization 
relaxes as if the initial state was the completly ordered 
state and one has 



m(n,t) = - ^Ji(4i) 



(16) 



for n > t. The local magnetization reaches the constant 
value —1/2 with corrections of order so that the 

derivative <P (v) is essentially vanishing for large n. This 
is exactly what is seen from the asymptotic behavior of 
the Bessel functions and their derivatives, vanishing as 
exp[-A(i;)n] with \{v) > 0. 

Inside the light-cone {v < 1), with the help of the 
asymptotics for 1 [0 



¥I7tan?C0sV' 

SHlMsinV 



(17) 



where tp = i^(tan /3 — /3) — 7r/4, one obtains for the deriva- 
tive <(t;): 



^ri(«) = \/l-V^ = <P (v) 



(18) 



which is n independent due to the exact cancellation of the 
sin and cos terms in (p^. Finally, by simple integration 
we obtain m(n, t) — <P{n/t) with the scaling function 



1/2 



V < -1 



^('") = < -i[w\/l -v^ + arcsinu] -1 < 1; < 1 (19) 
-1/2 V > 1 

This has to be compared with the XX chain result 
— ;|arcsinu and ±1 outside the causal region. Contrary 
to the XX chain which has a conserved dynamics (the 
total z-component of the magnetization is a constant of 
motion) in the Ising quantum chain there is a transient 
regime where the local magnetization relaxes faster than 
toward the stationary value ±1/2 and then only the 
residual kink spreads as in the XX chain. The results 
are shown in figure 1 where the inset describes the initial 
transient regime. 



A - 




Fig. 1. Nonequilibrium transverse magnetization scaling 
function for the Ising quantum chain. The analytical expres- 
sion (^9]) and the numerical results are indistinguishable. In 
the inset, the transient regime is shown for times smaller than 
t = 2. The magnetization relaxes toward the value ±1/2. 



4 Droplet-like initial state 

Let us consider now the following initial state 



TTTD-TTT 



(20) 



with I 4) = I i)®"^, that is a droplet of L down spins inside 
a bath of up spins, both interacting and evolving with the 
quantum Hamiltonian H. This can be considered as a toy 
model for a quantum system (the midle part) coupled to 
some environnement (the external part), both governed by 
the same microscopic interactions and one can study how 
the system part relaxes due to the coupling to the external 
degrees of freedom. We start with a one spin droplet within 
the Ising model. In this case, using ( p^ ) together with (pO[), 
the Ising transverse magnetization is given by 



m{l, t) 



At 



(21) 



where ^;(t') is the function introduced in the previous sec- 
tion. For t; > 1 the magnetization is dominated by the first 
two terms since then the function (p^ (I / 1) is exponentially 
small. On the other hand for u < 1, after the faster re- 
laxation toward 1/2 a scaling regime emerges for the local 
excess magnetization, that is 



m^il^t) = m{l,t) - 1/2 = t-'^<P (-) 



(22) 



where <P (v) is given by (18). 

More generaly, for a droplet of size L, the transverse 
magnetization at time t is given by 



L/2 



J J2 1>lMl-k)/t) 



(23) 



k=~L/2 
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In the light-cone the excess magnetization is given to the 
dominant order in by 

/L/2 A/t+L/2t 
<p'{{l~k)/t)dk= / <p'{u)du 
-L/2 Jl/t~L/2t 

(24) 

so that finally in the scaling regime I L we have simply 



2L 



(25) 



For / 0{L), with J„(z) y/^/nz cos(z-n7r/4-7r/4) for 
large the local magnetization excess is simply —2L/TTt 
plus subdominant corrections. Then the total magnetiza- 
tion remaining inside the initial droplet region is 



2£_ 

nt 



(26) 



which is exactly what was obtained in ref. ||] for the XX 
quantum chain. As in the XX case, where the total mag- 
netization is conserved, something similar happens in the 
Ising case. If one considers the total magnetization excess 
at time t, given by the integral over the whole space, we 
simply have a constant: 



m''{l,t)dl -Li <P {v)dv 



(27) 



where — L is the residual excess magnetization after the 
initial transient regime since the up(down) domain relaxes 
locally toward l/2(— 1/2). This means that after the ini- 
tial loss of magnetization, which takes place on micro- 
scopic time scales of order t ~ 1/h = 1, the dynamic 
is conservative. We have a conservative deviation around 
the stationary value. In fact, one can show that the local 
magnetization excess satisfy a lattice continuity equation 
dtm'^{l,t)+j{l,t) — j{l — l,t) — 0, with the current density 
j{l,t) given by 

L/2 

jil,t) = -2 J2(/-fc)-i(2<)J2(/-fe)(2i) 

k=-L/2 

-•^2(;-fe)-2(2i) J2(/-/c)+i(2i) (28) 

which is related to the expectation value of Ff^^F^ oc 
ff_]^(Tf . In the continuum limit, the current is simply ex- 
pressed as j{x,t) = jm'^{x,t). 

Although, the total residual magnetization in the sys- 
tem part was calculated in ref. |^ for the XX chain, the 
scaliiig profile was not considered. Using equations ( pT|) 
and (|l^ ) we have for all values of the transverse field h 



m^ilt) = m{l,t) 



L/2 

l = -2 ^ Jl_,{t). (29) 

fc=-L/2 



In the light-cone, together with the asymptotic expressions 
for the Bessel functions, we obtain in the scaling regime 
t> i:s> F 

m^{l,t)^-'F(-] (30) 



t 



with the scaling function 
W{v) = 



(31) 



One may verify that the integral over the whole space of 
the local magnetization gives back —2F as it should be 
for the conserved dynamics system under consideration. 
The scaling functions for both Ising and XX chains are 
presented in figure 2, where the analytical results are com- 
pared with numerics. 
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Fig. 2. Scaling functions for the critical Ising and XX{mset) 
quantum chains. The oscillations in the XX case are finite size 
effects due to the initial droplet size. One can see that as the 
time is increased the numerical results are closer and closer to 
the analytical scaling function. This can be seen more evidently 
at the boundaries v = ±1 for the Ising chain. 



5 General initial z-state 

Clearly, for translation invariant Hamiltonians, equa- 
tion (|l^) giving the transverse magnetization is a discret 
convolution product: 



1(1, t)= J2 Ftil-k)S{k)^{Ft*S)il) 



(32) 



/e— — oo 



with S{k) = (<P'|a^|if'). The kernel Ft{l) is given in the 
continuum limit by Ft{l) = ^/(f) with 



v\ < 1 
v\ > 1 



(33) 



where the k — 1 refers to the Ising case and k = to the 
XX chain. The local magnetization m(x,t) = mt{v), is 
then expressed as the convolution product 



mtiv) ^ (St * f)iv) 



(34) 
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with St{v) = S{tv). For the kink hke initial state, we have 
St{v) — —sgn{v) — l — 2H{v), where H{v) is the Heaviside 
function. For the droplet-like initial state, 

where n{x) is the caracteristic function of the interval 
[—1/2, 1/2]. In the long time limit, i ^ L, we have St{v) — 
1 — ^^S{v), so that we recover very simply the results of 
the previous section. 

For a general initial state in the z-direction, the Fourier 
transform of equation ( ^2[) is 

mt{q) ^ St{q)f{q) (36) 
with ^ 

fniq) = - [ (1 - v^r^^/^e-^'^'^^dv (37) 

71" J-1 

SO that the kernels in Fourier space are simply given by 

/o(g) = M27rq) (38) 
for the XX chain and 

fi{q) = ^^^^ = ^iM^^q) + M2nq)) (39) 

for the Ising chain. By inverse Fourier transform, it is pos- 
sible to obtain the desired magnetization profile in direct 
space. For exemple, if we consider the modulated initial 
state S{x) — cos{2ttx/L), with L ^ 1, we obtain in the 
long time regime t ^ L 

m{x,t) = cos{2'!Tx/L)f^{2TTt/L) , (40) 

that is the modulation does not spread with time but only 
the amplitude is decreasing as t~^/^~'^. 

Finally, for a homogeneous initial state with S{x) = 
m(0), it is easy to see that m(t) — m(0) in the XX case, 
since the dynamic is conservative, and m{t) = l/2r7i(0) 
for the Ising chain. 

6 Summary 

We have calculated for the critical Ising and XX quan- 
tum chains the nonequilibrium transverse magnetization 
profiles for kink- like and droplet-like initial states. In both 
cases at large times the magnetization profiles exhibit scal- 
ing forms t^^F{l/t) which have been obtained analytically 
and are in excellent agreement with the numerics. The 
two systems show essentially the same features even if the 
dynamics of the transverse magnetization are very differ- 
ent, conservative for the XX chain and nonconservative 
for the Ising model. In the Ising case there is a transient 
regime where the initial magnetization relaxes toward the 
homogeneously initialy magnetized state stationary value 
±1/2. After this initial regime, the system evolves as if 
the dynamics of the residual transverse magnetization was 
conservative. The long time relaxation of the transverse 
magnetization starting with a general initial z-state is ex- 
pressed very simply in terms of a convolution product of 
the initial distribution with a response kernel obtained 
analytically for both XX and Ising chain. 



Appendix: Time evolution 

The diagonalisation of the Hamiltonian (||) leads to 

^ = *E|«- (41) 

In the Fermi operator representation, with rjq = 1/2(7^ -|- 
17^) and ryj the hermitian conjugate, one obtains H — 
Tlq ^qVlVq~{^/'^) J2q ^q- The time evolution of the Chfford 
operators is given by U^{t)'-fqUq{t) with 

Uq{t) = exp (^7^,^) - cos ^ + ^^q sin ^ , (42) 

which leads to equation (^). Since {7^,7^/} = 2SijSqq', 
we can write equivalently for equation (^ 7g(^) — 
J2'j=i{lq\lqit))lqj whcrc the symbol (.|.) means the half 
of the anticommutator. 

The time evolution of the Fs is then expressed as 

A' (^ ) - I] 09 (fc) cosieqt)-fl +<j,q{k) sixv{eqt)-il (43) 
q 

and 

A'(i) - E -^«(^) sin(egi)7g + ^q{k) cos{eqt)-il (44) 
q 

with initial values rl{0) = Y.q<l^q{k)ll and ^^(O) = 
J2q''J-'qi^)lq- Reinjecting in this expressions the inverse 
transforms 7^ = Y.k^q^^)^k and 7^ = Y.k^qi^)rk one 
finally arrives at equations (^ with components (|8[). 

Formally, since the anticommutators {r'^,/^} = 
25ij5kj are all proportional to the identity operator, the 
set {-Tfc} forms an orthonormal basis of a 2L-dimensional 
linear vector space £ with inner product defined by (.|.) = 
^{., .}. Hence, every vector X & £ has a unique expansion 
X = Y.^,k{^i\X)^l. The string expression XiX2...X„, 
with Xj g £, is a direct product vector of the space 
£1 ® £2 ® ■■■ ® £n which decomposition is 

M...x„= {ri}^\x{)...{rii\x^)rii..Tii 

ii .ki , . . .,in,kn 

(45) 

Using this formalism, the local magnetization at time t is 
given by 

= E ini\rnt)){ni\r,\t))riirii{46) 

ii,fci,i2,fe2 

which is our starting point. One has then to con- 
sider the simple time-independant expectation values 
(iZ'lr'^jr'^^ jiZ'), which are easilly obtained in the spin basis 
using the Jordan- Wigner expressions 
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